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Abstract 

Analysis  of  Argo  float  trajectories  at  1000  m  and  temperature  at  950  m  in  the  North  Atlantic  between  November  2003  and  January 
2005  demonstrates  the  existence  of  two  different  circulation  modes  with  fast  transition  between  them.  Each  mode  has  a  pair  of  cy¬ 
clonic  -  anticyclonic  gyres.  The  difference  is  the  location  of  the  cyclonic  gyre.  The  cyclonic  gyre  stretches  from  southeast  to  north¬ 
west  in  the  first  mode  and  from  the  southwest  to  the  northeast  in  the  second  mode.  The  observed  modes  strongly  affect  the  heat  and 
salt  transport  in  the  North  Atlantic.  In  particular,  the  second  mode  slows  down  the  westward  transport  of  the  warm  and  saline  water 
from  the  Mediterranean  Sea. 
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1  Introduction 

Although  the  North  Atlantic  is  a  subject  of  in¬ 
tensive  theoretical  and  observational  studies  for  many 
years ,  knowledge  on  the  large  -  scale  mid  -  depth 
circulation  and  especially  on  its  variability  is  quite 
limited.  There  are  several  attempts  to  construct  the 
large  -  scale  mid  -  depth  circulation  pattern  for  the 
whole  North  Atlantic  using  climate  data  (  Schmitz 
and  McCartney,  1993;  Reid,  1994;  Lozier  et  al.  , 
1995 )  and  for  regional  seas  in  the  North  Atlantic 
with  higher  resolution  using  subsurface  float  (  SO- 
FAR,  RAFOS,  ALACE,  PALACE,  SOLO)  trajec¬ 
tories  directly  or  combining  with  hydrological  obser¬ 


vations  (  Lavender  et  al.  ,  2000 ;  Lavender  et  al.  , 
2005;  Zhang  et  al.  ,  2001;  Bower  et  al.  ,  2002). 
The-basin  scale  and  subbasin-scale  mid-depth  circu¬ 
lation  and  its  variability  can  also  be  potentially  ex¬ 
tracted  from  the  Argo  float  data  (  Argo  Science 
Team  ,  2001 ) .  However,  analysis  of  Argo  data  is  not 
a  trivial  problem  because  the  observation  is  sparse 
(  up  to  800  km  spatial  gap )  and  noisy. 

2  Data  and  quality  control 

Between  November  2003  and  January  2005 , 
over  56  000  float  days  (  cumulative )  of  data  were 
collected  in  the  North  Atlantic  ( 10°  ~60°N)  in  gen¬ 
eral  at  three  parking  depths ;  1  000 ,  1  500  and 
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2  000  m.  The  floats  parking  at  2  000  m,  depths 
shallower  than  1  000  m ,  and  unknown  depths  are  ex¬ 
cluded  from  the  analysis.  Temperature  at  950  m  and 
trajectories  at  1  000  and  1  500  m  are  extracted  from 
all  the  existing  Argo  floats.  The  data  from  1  000  and 
1  500  m  were  grouped  together  to  represent  the  mid¬ 
depth. 

The  measurement  cycle  of  an  Argo  profiling 
float  includes  four  stages;  ascending,  surface  drif¬ 
ting,  diving  and  deep  drifting.  The  Argo  float  can 
only  get  its  position  fixings  while  it  ascends  to  the 
sea  surface.  The  vector  between  two  consecutive  sur¬ 
face  positions  during  the  deep  drifting  divided  by  the 
time  interval  is  taken  as  the  mid  -  depth  velocity 
vector.  When  the  Argo  float  is  diving,  ascending 
and  drifting  below  the  sea  surface ,  no  data  can  be 
transmitted  to  the  ground  stations  in  real  time.  Ve¬ 
locity  field  after  the  first  step  analysis  shows  noisy 
circulation  patterns  (  Figs  1  a ,  and  b )  with  large  spa¬ 
tial  gaps  (  from  230  to  800  km ) . 

Uncertainty  in  the  Argo  float  data  causes  errors 
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in  the  velocity  field.  First ,  the  data  extracted  from 
the  floats  parking  at  two  different  levels ;  1  000  and 
1  500  m  and  grouped  together  to  represent  the  mid  - 
depth  (1000  m).  This  neglects  the  vertical  shear. 
Second ,  the  vertical  shear  causes  increase  or  de¬ 
crease  of  the  distance  between  the  points  of  ascend¬ 
ing  from  and  diving  to  the  parking  depth.  Thirdly , 
the  sequence  of  float  trajectory  segments  only  approx¬ 
imates  the  real  Lagrangian  paths.  Fourthly,  prelimi¬ 
nary  computations  (  not  included  here  to  be  pub¬ 
lished  in  a  separate  paper)  show  that  high  resolution 
elements  of  circulation  in  the  western  North  Atlan¬ 
tic  ,  such  as  the  northern  recirculation  gyre  and  the 
Deep  Western  Boundary  Current  (DWBC)  are  also 
revealed  by  the  Argo  floats.  For  example,  Figs  la 
and  b  clearly  show  the  existence  of  DWBC.  Howev¬ 
er,  such  a  resolution  is  not  available  for  the  whole 
North  Atlantic.  The  high  -  energetic  mesoscale  ed¬ 
dies  and  narrow  boundary  currents  are  classified  as 
“  noise  ”  and  removed  from  the  analysis.  Fifthly , 
there  are  large  spatial  gaps  in  Argo  float  trajectories. 
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Fig.  1.  Circulation  velocities  (tiny  arrows)  estimated  from  the  original  Argo  float  tracks  at  1000  m  for  December  2003 — 
Mar  2004  (a)  and  August  2004 — November  2004  (b).  The  figure  scale  is  given  for  tiny  arrows.  Red  arrows  are  circulation 
velocities  obtained  by  averaged  over  appropriate  bins  ( red  lines ) .  The  magnitudes  of  these  arrows  are  multiplied  by  the  scale 
equaled  20  for  better  visualization  of  schematic  circulation  patterns.  i\  ,  /  2  and  /”3  are  boundaries  of  the  computation  area. 
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The  temperature  observation  has  higher  quality 
than  the  velocity  observation  since  the  former  has 
( 1  )  higher  resolution  and  ( 2  )  less  spatial  gaps. 
Figure  2  shows  typical  monthly  coverage  of  observa¬ 


tions.  Interested  readers  can  find  the  detailed  dis¬ 
cussion  on  the  navigation  errors  and  measurement  er¬ 
rors  caused  by  temperature  sensors  in  http;//argo. 
jcommops.  org. 
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Fig.  2.  Monthly  temperature  at  950  m  for  Feburary  2004  (a)  and  December  2004  (b).  Circles  are  points  where 
temperature  was  measured.  The  zero  flux  boundary  condition  is  used  at  _T\  ,  r2  and  r3  to  reconstruct  temperature  pat¬ 
terns  (see,  explanation  in  Chu  et  al.  ,  2003a). 


Three  quality  control  steps  are  performed  to  i- 
dentify  and  remove  temperature  profiles  corrupted  by 
large  measurement  errors.  (  1  )  Explicitly  absurd 
profiles  were  removed  after  a  visual  inspection.  (2) 
A  portion  of  temperature  profiles ,  which  were  outside 
of  the  prescribed  accuracy  of  the  climatic  data  ( Lev- 
itus  et  al.  ,  1998)  ,  were  also  excluded  from  the  fur¬ 
ther  analysis.  ( 3  )  Temperature  snapshots  computed 
by  the  optimal  spectral  decomposition  (  OSD )  meth¬ 
od  (Chu  et  al.  2003a,  2004)  should  not  show  ex¬ 
plicit  outbreaks  in  temperature  structure.  Any  tem¬ 
perature  profile  contributing  to  such  outbreaks  was  a 
subject  to  removing. 

A  parking  depth  for  each  Argo  float  was  extrac¬ 


ted  from  a  “meta”  file  as  the  variable  “PARKING_ 
PRESSURE”  (  http ;/ / www.  Argo.  ucsd.  edu  )  .  To 
control  the  parking  depth  the  variable  “ PRES ”  (a 
pressure  measured  along  the  float  trajectory)  from  a 
file  containing  float  trajectory  data,  was  used.  Most 
floats  launched  in  the  area  of  interest  have  measured 
this  variable. 

3  Methodology 

Two-step  analysis  is  used  in  this  study.  First , 
the  bin  technology  (Davis,  1991 )  is  used  to  get  the 
first-guess  mid-depth  velocity  (  temperature  )  field 
from  the  Argo  float  trajectories  (  temperature  )  at 
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1000  m  (950  m)  (see  Figs  la  and  b).  Second,  the 
optimal  spectral  decomposition  (  OSD )  method  (  Er- 
emeev  et  al.  ,  1992;  Ivanov  et  al.  ,  2001;  Chu  et 
al.  ,  2003a,  b,  2004;  Chu  et  al.  ,  2007)  is  used  to 
reconstruct  circulation  and  temperature  patterns  from 
the  first-guess  fields  with  higher  resolution. 

Any  field  (  temperature ,  salinity ,  or  velocity ) 
can  be  decomposed  into  generalized  Fourier  series  u- 
sing  the  OSD  method.  The  three-dimensional  field  is 
then  represented  by  linear  combination  of  the  prod¬ 
ucts  of  basis  functions  ( or  called  modes )  and  corre¬ 
sponding  Fourier  coefficients  at  point  x  of  the  inter¬ 
ested  area  (W, 

L 

T(x,t)  =T0+^Dl(t)Ul(x,t),  (1) 

1  =  1 

where  \  S\\  is  the  basis  functions ;  j  D:  (  t )  )  is  the 
spectral  coefficients ;  and  T0  is  the  mean  temperature 
over  the  interested  area  ( Ivanov  et.  al.  ,  2001 ;  Chu 
et  al.  ,2004).  The  least  square  method  is  used  to 
determine  the  spectral  coefficient  and  the  mean  tem- 
peralure  (  Thiebaux  and  Pedder ,  1987 ) . 

If  a  rectangular  closed  ocean  basin  is  consid¬ 
ered,  the  basis  functions  are  sinusoidal  functions.  If 
a  realistic  ocean  basin  is  considered ,  the  basis  func¬ 
tions  are  the  eigenvalues  of  the  three  -  dimensional 
Laplace  operator  with  real  topography.  One  major 
benefit  of  using  the  OSD  method  is  that  the  boundary 
conditions  for  the  ocean  variables  ( temperature ,  sa¬ 
linity  ,  velocity )  are  always  satisfied. 

After  the  decomposition ,  the  three  -  dimension¬ 
al  field  is  represented  by  a  set  of  Fourier  coeffi¬ 
cients.  This  method  has  three  components  ;  (  1 )  de¬ 
termination  of  the  basis  functions ;  ( 2 )  optimal  mode 
truncation;  and  (3)  determination  of  the  Fourier  co¬ 
efficients.  Determination  of  basis  functions  is  to 
solve  the  eigenvalue  problem.  Chu  et  al.  (2003a, 
b)  developed  a  theory  to  obtain  the  basis  functions 
with  open  boundaries.  The  basis  functions  are  only 
dependent  on  the  geometry  of  the  ocean  basin ,  not 
dependent  on  the  oceanic  variables.  This  is  to  say, 


no  matter  which  variable  ( temperature ,  salinity ,  or 
velocity )  is  concerned ,  the  basis  functions  are  the 
same ,  and  can  be  predetermined  before  the  data  a- 
nalysis.  For  data  without  error,  the  more  the  modes, 
the  more  the  accuracy  of  the  processed  field.  For  da¬ 
ta  with  error,  this  rule  of  the  thumb  is  no  longer 
true.  Inclusion  of  high  -  order  modes  leads  to  in¬ 
creasing  error.  The  Vapnik  variational  principal 
(Vapnik,  1982)  is  used  to  determine  the  optimal 
mode  truncation.  After  the  mode  truncation ,  optimal 
field  estimation  is  to  solve  a  set  of  a  linear  algebraic 
equation  of  the  Fourier  coefficients.  This  algebraic 
equation  is  usually  ill  -  posed.  The  rotation  method 
(Chu  et  al.  ,  2004)  is  developed  to  change  the  ma¬ 
trix  of  the  algebraic  equation  from  ill  -  posed  to 
well  -  posed  such  that  a  realistic  set  of  the  Fourier 
coefficients  are  obtained. 

Thus ,  the  OSD  method  has  two  major  features. 
First,  a  special  constraint  (Vapnik,  1982)  is  used 
to  determine  the  optimal  truncation  for  to  compromise 
between  spectral  representation  and  ability  to  fit  the 
observations.  For  number  of  basis  functions  more 
than  the  optimal  truncation ,  the  model  contains  too 
much  contributions  from  noise.  For  number  of  basis 
functions  less  than  the  optimal  truncation ,  the  model 
resolution  is  too  coarse  to  fit  the  observations. 

The  optimal  truncation  determines  the  real  spa¬ 
tial  resolution  of  the  reconstructed  field  and  depends 
on  the  noise  level  and  degree  of  inhomogeneity  of  ob¬ 
servation  coverage.  It  fills  large  spatial  gaps  in  ob¬ 
servation  coverage.  Interested  readers  are  referred  to 
Chu  et  al.  (2003a,  b,  2004;  Chu  et  al.  ,  2007) 
where  all  the  technical  details  and  features  of  the 
OSD  method  are  explained.  After  the  spectral  coeffi¬ 
cients  and  mean  teroperature  are  determined  from  the 
observational  data  T(xp)  ,  the  temperature  field  with 
regular  temporal  -  spatial  grids  is  obtained  using  Eq. 
(1). 

A  vector  field  (  such  as  velocity )  at  the  point  x 
of  the  interested  area  is  represented  by 
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N 

U(x)  =  CV%(x)  +  x  [kVn(x)]  + 

n  =  1 
M 

-0m(x),  (2) 

m  =  1 

where  V/ ( x )  and  <Pm(x)  are  the  basis  functions;  An 

and  Bm  are  the  spectral  coefficients ;  k  is  the  unit  vec¬ 
tor  in  the  vertical  ( upward  positive )  direction ;  V  is 
the  horizontal  gradient  operator ;  V/  (  x  )  and  C  are  the 
harmonic  function  and  its  coefficient  for  the  multicon- 
nected  area  of  the  interested  area  with  three  bounda¬ 
ries  l\  ,  r2,  r3  (  Fig.  la).  All  other  notations  are 
the  same  as  for  temperature.  The  optimal  mode  trun¬ 
cations  Mopt  and  /Vorit  are  determined  using  Vapnik’s 
( 1982 )  technology.  After  the  basis  functions  [  V/ 
(x)  ,  Om(x)  I  re  =  1 ,  ••• ,  N opt ;  in  =  l,---,  Mopt  ]  the 
harmonil  function  and  its  coefficient  and  V/  ( x  )  , 
C  ]  are  determined  from  the  observational  data 

U(x)p  =  [  U\  (  xp  )  ,  U2  ( xp  )  ]  collected  within  some 
time  period ,  the  velocity  vector  field  with  regular 
temporal  -  spatial  grids  is  obtained  using  Eq.  (2). 

Comparison  among  OSD  ,  optimal  interpolation , 
and  smoothing  -  spline  (  Chu  et  al.  ,  2003a,  b) 
shows  great  advantage  of  the  OSD  method  in  filling 
large  gaps  (  spatial  and  temporal )  in  observational 
coverage  and  removing  noises  with  various  scales  ef¬ 
fectively  from  the  data.  These  features  are  the  key  to 
select  the  OSD  method  for  the  analysis  of  the  Argo 
data. 

4  Estimation  of  positioning  errors 

As  we  mentioned  previously,  no  data  can  be 
transmitted  to  the  ground  stations  in  real  time  below 
the  sea  surface.  However,  the  vertical  geostrophic 
shear  may  be  strong  especially  in  the  west  boundary 
current.  Therefore,  estimation  of  positioning  errors 
is  needed  in  calculating  the  intermediate  -  depth  ve¬ 
locity  from  Argo  trajectory  data.  Park  et  al.  (2005  ) 
have  assumed  that  a  surface  Argo  float  trajectory 
consists  of  large  -  scale  drift ,  inertial  motion ,  and 


measurement  noise.  For  temporal  scales  longer  than 
two  months ,  the  inertial  motions  may  be  considered 
as  spatially  uncorrelated  “noise”  and  effectively  re¬ 
moved  by  the  OSD. 

Systematic  float  shifts  caused  by  large  -  scale 
currents  in  upper  ocean  layer  are  estimated  following 
Chu  et  al.  (2007).  A  location  of  the  pth  float  x  at 
the  parking  depth  H  0  is  specified  by  deduction  of  the 
“effective”  displacement  from  a  location  of  the  same 
float  xp  after  ascending  or  before  descending  to  the 
parking  depth ; 


where  St  =  t  —  t'  ,  is  the  ascending/ descending 
time ;  w  and  is  a  velocity  of  ascending  or  descend¬ 
ing.  A  phenomenological  model  is  used  for  U  ; 

v  =  | Up  exp (  -  crpz )  ,  z^Hp,  (4) 
'  z>Hp, 

where  Usp  and  0  are  the  surface  and  mid-depth  ve¬ 
locities  estimated  from  the  surface  and  subsurface 
tracks  of  the  pth  float ,  respectively ;  //  is  the  low 
boundary  of  the  oceanic  thermohaline  estimated  from 
temperature  profile  of  the  same  float ;  crp  =  -  1  / 

(2 Hp)  log  (U2/U2).  For  Argo  floats  launched  in 
the  area  of  interest  we  have 

UsS  «  15  -  20  cm/s,  U  «  5  -  10  cm/s, 

H  «  300  m,  St  =  10  h.  (5) 

Substituting  values  of  Eq.  (  5  )  into  Eqs.  (  3  ) 
and  ( 4 )  yields  that  on  an  average  the  position  error 
for  the  Argo  float  launched  in  the  area  of  interest 
should  be  less  than  6  km.  That  results  in  the  system¬ 
atic  drift  for  floats  less  than  0.  6  cm/s.  Therefore, 
although  the  position  error  of  each  float  was  corrected 
by  the  method  described  above ,  the  upper  ocean 
drifting  effects  does  not  contribute  too  much  into  the 
reconstruction  error. 

5  Results 

The  two  different  circulation  modes  (  Figs  1  a 
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and  b)  are  barely  recognized  from  the  Argo  float  ve¬ 
locity  using  the  bin  technology.  The  arrows  are  the 
mean  velocities  averaged  over  appropriate  bines 
shown  in  the  same  figures.  The  circulation  patterns, 
robust  to  variations  of  bin  sizes ,  show  well  -  known 
circulation  features  identified  from  early  analyses  on 
the  RAFOS  float  trajectories  in  separate  regions  of 
the  North  Atlantic.  For  example,  across  the  Mid-At¬ 
lantic  Ridge  transport  of  cooler  and  fresher  Labrador 
Sea  water  at  about  50"  N  (  Paillet  et  al.  ,  1998  )  is 
clearly  seen  in  Fig.  la.  Note  that  the  bin  technology 
is  unable  to  remove  stochastic  noises  from  the  obser¬ 
vations.  Davis  (  1991  )  also  pointed  out  numerical 
biases  of  this  technology.  Therefore ,  the  results  from 
the  bin  technology  (  Figs  1  a  and  b )  are  crude  as¬ 
sessments  of  the  reality. 

Different  temperature  patterns  should  corre¬ 
spond  to  these  circulation  modes  since  they  are  very 
distinct  features.  To  confirm  this ,  we  reconstructed 
14  monthly  temperature  fields  between  December 
2003  and  January  2005  with  the  optimal  truncation  of 
38 ,  temperalure  of  6.  50  °C  and  the  basis  function  on 
1°  X  1°  grid.  These  computations  uniquely  demon¬ 
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strate  the  existence  of  two  different  temperature  pat¬ 
terns  for  February  2004  (Fig.  2a)  and  December 
2004  (Fig.  2b). 

The  principle  distinction  between  these  patterns 
appears  at  least  in  cross-Atlantic  transport  of  heat 
from  the  Mediterranean  Sea.  The  first  mode  implies 
the  existence  of  advection  pathways  across  the  Mid  - 
Atlantic  Ridge ,  which  may  accelerate  westward  mo¬ 
tion  of  the  Mediterranean  water  in  the  North  Atlan¬ 
tic.  The  second  mode  slows  down  westward  propaga¬ 
tion  of  warm  and  salt  water  from  the  Mediterranean 
Sea.  Comparison  of  monthly  temperature  patterns 
consecutively  from  September  to  December  2004 
( the  second  mode  period ,  not  shown  here  )  clearly 
demonstrates  a  largest  portion  of  Mediterranean  water 
slowly  moves  along  the  Mid-Atlantic  Ridge  at  speed 
of  about  1.0  ~  1.5  cm/ s.  If  the  mid-depth  tempera¬ 
ture  is  taken  as  a  tracer,  the  velocity  can  be  roughly 
estimated  from  two  consecutive  temperature  fields  on 
the  base  of  displacements  of  6 ,  7 ,  8 ,  9  and  10  T! 
isotherms  (  Figs  3  a  and  b )  where  the  arrows  show 
the  mean  direction  of  isotherm  displacements  during 
three  month  periods. 
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Fig.  3.  The  cartoons  of  mid-depth  circulation  drawn  from  the  analysis  of  isotherm  displacements  for  December  2003 — March 
2003  (a)  and  September  2004 — December  2004  (b).  Dashed  and  solid  lines  are  monthly  isotherms  for  Dectember  2003- — 
March  2004  in  a  and  for  September  2004 — December  2004  in  b.  The  gray  arrows  schematically  show  water  pathways. 
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Spatial  distributions  of  both  modes  are  recon¬ 
structed  for  7Vopt  =  24 ,  Mop  =  3,  C  =  3.  2  x  107  and 
the  basis  functions  calculated  on  1°  x  1°  grid.  Two 
distinct  circulation  modes  and  a  transition  between 
them  are  identified  from  these  computations 


(  Figs  4  a ,  b  and  c ) .  The  first  mode  represents  a 
pair  of  cyclonic  -  anticyclonic  gyres  and  two  gyres  of 
smaller  size  in  the  Labrador  Sea  and  the  eastern 
North  Atlantic  (  Fig.  4a ) .  This  mode  dominates 
from  December  2003  to  March  2004. 
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Fig.  4.  Mean  stream-function  at  1000  m  for  December  2003 — March  2004  (a)  April  2004 — July  2004  (b)  and  August 
2004  November  2004  ( c )  computed  by  the  OSD  method.  Contours  are  labeled  in  units  of  107  cm2/ s.  The  kinematic 
boundary  and  open  boundary  conditions  used  at  _T\  and  T2 ,  and  P3 ,  respectively,  are  explained  in  Chu  et  al.  (2003). 
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The  anticyclonic  gyre  has  an  intensive  core  in 
the  western  part  of  the  North  Atlantic  with  two  bran¬ 
ches  stretched  to  the  northeastern  and  southeastern 
Atlantic  (  Fig.  4a ) .  The  center  of  the  core  is  not 
fixed  in  the  same  location.  It  may  shift  to  the  north¬ 
east  by  up  to  1000  km.  The  maximum  velocity  is 
about  4  ~  5  cm/ s  in  the  core  and  0.  5  ~  1  cm/ s  in  the 
branches  of  the  gyre.  The  cyclonic  gyre  stretches 
from  the  southeast  to  northwest  up  to  25  °N  with  a  ve¬ 
locity  about  1  ~  2  cm/ s. 

The  second  mode  represents  a  pair  of  cyclonic 
and  anticyclonic  gyres  on  basin  scales  (  Fig.  4c ) . 
Flowever,  in  this  case  the  cyclonic  gyre  stretches 
from  the  southwest  to  the  northeast  up  to  43  °N  and 
the  boundary  between  these  gyres  lies  along  the 
Mid  -  Atlantic  Ridge.  The  velocity  reaches  3  ~ 
4  cm/ s  in  the  anticyclonic  gyre  and  2  ~  3  cm/ s  in 
the  cyclonic  gyre. 

Reid  ( 1994 )  identified  the  second  mode  as  the 
steady  state  of  the  North  Atlantic  due  to  its  high  re¬ 
currence  in  data.  Our  computations  indicate  that  this 
mode  is  one  of  the  two  possible  North  Atlantic 
modes.  The  two  modes  correspond  to  opposite  mean 
values  of  the  North  Atlantic  oscillation  (  NAO )  index 
(www.  cru.  uea.  ac.  uk)  ;  negative  value  (  -  0.  81 ) 
for  the  first  mode  (  December  2003  to  March  2004 ) 
and  positive  value  ( 0.  27 )  for  the  second  (  August — 
November  2004  )  .  Explicit  and  fast  transitional 
inter  -  mode  is  clearly  seen  between  April  2004  and 
July  2004  (Fig.  4b). 

It  is  noted  that  a  cold  center  is  identified  near 
20°N,  55°W  in  February  2004  (Fig.  2a)  with  large 
spatial  gaps  in  Argo  data.  Comparing  with  the  recon¬ 
structed  stream-function  from  the  position  data  for 
December  2003 — March  2004  (Fig.  4a)  ,  a  cyclonic 
circulation  is  associated  with  the  cold  center.  This 
confirms  that  the  southwestern  cold  center  shown  in 
Fig.  2a  is  physical  and  not  caused  by  the  large  spa¬ 
tial  gaps  in  Argo  float  temperature  data.  Here ,  the 
four  months  "  positioning  data  (  December  2003 — 


March  2004  )  is  used  due  to  the  data  sparseness. 
This  also  reflects  the  usefulness  of  the  OSD  method. 

6  Sensitivity 

Our  estimations  demonstrate  that  neglecting  the 
vertical  shear  between  1000  and  1500  m  does  not 
cause  large  bias  of  the  reconstructed  circulation  pat¬ 
terns  (Figs  4  a,  b  and  c).  To  confirm  this,  we  re¬ 
constructed  circulation  at  1000  m  or  1500  m  sepa¬ 
rately  in  certain  regions  with  the  floats  parking  at  the 
same  depth,  and  compared  the  circulation  patterns 
in  that  part  of  the  North  Atlantic  where  these  areas 
intersect.  We  did  not  find  considerable  differences 
between  the  circulation  patterns  at  basin  -  scales. 

Mode  truncations  at  24  for  velocity  and  38  for 
temperature  are  robust  to  variations  of  the  length  of 
samples  ( number  of  observations )  and  the  level  of 
noise  distorting  the  Argo  observations.  The  Euclide¬ 
an  norms  ||  A  ||  ^,  and  ||  D  ||  ls  increase  only  for 

5%  ~  6%  and  1%  ~  3%  with  the  increase  of  ( N , 

L)  to  (50,  100)  ,  and  vary  little  even  if  approxi¬ 
mately  10%  ~  15%  of  the  data  are  excluded  from  the 
reconstruction  process.  Robustness  of  the  reconstruc¬ 
ted  fields  is  an  additional  evidence  for  the  recon¬ 
struction  reliability. 

7  Computation  error 

Velocities  computed  along  original  Argo  float 
trajectories  (  non  -  filtered  )  are  distorted  by  quasi  - 
isotropic  non  -  Gaussian  white  noise  with  zero  mean , 
standard  deviation  (  5.312  cm/s  ),  skewness 

(0.  034)  ,  and  kurtosis  (6.  34)  ,  mixed  with  spatial¬ 

ly  correlated  (  red )  noise  with  zero  mean ,  standard 
deviation  ( 0.  049  cm/ s )  ,  and  the  correlation  radius 
(150  ~ 200  km).  The  white  noise  represents  the  re¬ 
sidual  between  the  original  observations  and  the 
same  observations  averaged  over  4"  x  4"  bins.  The 
red  noise  is  the  observational  velocity  minus  the  sum- 
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mation  of  the  reconstructed  velocities  and  white 
noise.  Temperature  observations  are  distorted  by 
quasi  -  Gaussian  white  noise  with  zero  mean  and 
standard  deviation  of  0.  1  °C . 

The  “noise”  in  the  velocity  observations  can  be 
effectively  filtrated  by  the  OSD  methad  using  the  ro¬ 
tation  method  (  Chu  et  al.  ,  2004)  if  the  number  of 
observations  is  greater  than  1800.  This  requirement 
can  be  reached  when  four  months  data  are  used. 
Similarly ,  the  "  noise "  in  the  temperature  field  can 
be  filtered  out  if  the  number  of  observations  is  grea¬ 
ter  than  400.  It  is  possible  to  use  the  data  within  only 
one  month  time  period. 

Our  estimations  also  show  that  the  OSD  methad 
reduces  the  noise  to  signal  ratio  for  the  Argo  data  to 
0.  25  ~  0.  30  (  velocity )  and  0.  07  ~  0.  1  ( tempera¬ 
ture  ) .  With  these  noise  levels ,  the  reconstructed 
circulation  and  temperature  fields  are  quite  realistic 
(Chu  et  al.  ,  2003;  Chu  et  al.  ,  2004). 

Errors  in  estimating  the  stream-function  and 
temperature  are  computed  using  the  technology  of 
“laminar  ensemble”  (  Turchin  et  al.  ,  1971).  The 
reconstruction  error  for  temperature  is  less  than 
5%  ~  7%  over  the  whole  North  Atlantic.  The  maxi¬ 
mum  standard  deviations  of  the  error  for  the  stream- 
function  is  up  to  1.  6  x  107  cmVs  in  eastern  part  of 
North  Atlantic  and  in  an  area  centered  at  intersection 
of  20  °N  parallel  and  the  Mid  -  Atlantic  Ridge.  Com¬ 
paring  this  error  with  values  of  the  reconstructed 
stream-function  of  1.  0  X  108  cmVs  for  the  first  area 
and  4.  0  X  107  cm  /s  for  the  second  area,  we  can 
conclude  that  the  reconstructed  circulation  patterns 
are  statistically  significant  and  difference  between 
the  modes  is  not  due  to  noisiness  and  sparseness  of 
the  data  and/ or  unsatisfactory  mathematical  analy¬ 
sis. 

8  Conclusions 

The  present  study  demonstrates  the  usefulness 


of  the  Argo  float  data  for  identification  of  basin  and 
subbasin-scale  mid-depth  circulation  in  the  North 
Atlantic.  The  velocity  field  is  quite  noisy  and  sparse 
when  the  bin  technology  is  used  (  see  Figs  1  a  and 
b).  Two  different  circulation  modes  and  fast  transi¬ 
tion  between  them  are  identified  using  the  optimal 
spectral  decomposition  method.  The  second  mode  was 
detected  by  Reid  ( 1994 ,  see  Fig.  8e  from  his  pa¬ 
per)  using  the  climate  data  arrays.  This  paper  dem¬ 
onstrates  that  the  North  Atlantic  circulation  can 
switch  between  two  modes  rather  than  “this  large  - 
scale  pattern  of  mid  -  depth  flows  does  not  vary  too 
widely”  (Reid,  1994).  The  transition  time  between 
the  modes  is  less  than  three  months  and  there  is  no 
seasonal  periodicity  in  temporal  mode  variability. 
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